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Abstract 

We discuss the evolution of linear perturbations in a quintessence model in which the 
scalar field is non-minimally coupled to cold dark matter. We consider the effects of this 
coupling on both cosmic microwave background temperature anisotropies and matter per- 
turbations. Due to the modification of the scale of cold dark matter as p c = p^a~ 3+ ^, we 
can shift the turnover in the matter power spectrum even without changing the present 
energy densities of matter and radiation. This can be used to constrain the strength of the 
coupling. We find that the phenomenology of this model is consistent with current observa- 
tions up to the coupling power n c < 0.01 while adopting the current parameters measured 
by WMAP. Upcoming cosmic microwave background observations continuing to focus on 
resolving the higher peaks may put strong constraints on the strength of the coupling. 



1 Introduction 



Analysis of the Hubble diagram of high redshift Type la supernovae (SNe la) has discovered 
that the expansion of the Universe is currently accelerating [1]. In addition, combining 
measurements of the acoustic peaks in the angular power spectrum of the cosmic microwave 
background (CMB) anisotropy which indicate the flatness of the Universe [2] and the matter 
power spectrum of large scale structure (LSS) which is inferred from galaxy redshift surveys 
like the Sloan Digital Sky Survey (SDSS) [3] and the 2-degree Field Galaxy Redshift Survey 
(2dFGRS) [4] has confirmed that a component with negative pressure (dark energy) should 
be added to the matter component to make up the critical density today. 

The cosmological constant and/or a quintessence field are the most commonly accepted 
candidates for dark energy. The latter is a dynamical scalar field leading to a time dependent 
equation of state, uj^. Also, this scalar field has a fluctuating, inhomogeneous component 
in order to conserve the equivalence principle corresponding to the response of the new 
component to the inhomogeneities in the surrounding cosmological fluid [5]. 

Several new observational effects produced by the existence of quintessence are imprinted 
in the CMB anisotropies and the matter power spectrum when compared to models with 
the cosmological constant. The locations of the acoustic peaks in the CMB angular power 
spectrum are shifted due to their dependence on the amount of dark energy today and at 
last scattering as well as [6, 7]. Usually the Universe is dominated by the quintessence 
at late times (z ~ 0(1)), when the gravitational potential associated with the density 
perturbations is changed due to a time dependent u^. This enhances the CMB anisotropies 
at large angular scales by the integrated Sachs- Wolfe effect (ISW) [8] . Thus the amplitudes 
of both CMB angular and matter power spectra decrease at large scale for the fixed Cosmic 
Background Explorer Satellite (COBE) normalization compared to those in the cold dark 
matter model with a cosmological constant (ACDM) [9]. 

The possibility that a scalar field at early cosmological times follows an attractor-type 
solution and tracks the evolution of the visible matter-energy density has been explored 
[10]. This may help alleviate the severe fine-tuning associated with the cosmological con- 
stant problem. However, this still cannot explain the reason why the dark energy and the 
dark matter have comparable energy densities at present. Recently models considering the 
coupling of quintessence to dark matter have been investigated as a possible solution for 
this late time coincidence problem [11]. However, several of these models using a simple 
coupling can be ruled out by observational constraints [12]. 

These non-minimally coupled quintessence models have several different observational 
effects compared to the minimally coupled models. One of the most important effects 
is a different scaling of the cold dark matter (CDM) compared to that of CDM of the 
minimally coupled case. Since CDM scales as p c = p^a~ 3+ ^ where £ < 0, there will be 
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more CDM energy density at early epoch compared to the case with £ = (i.e. n c — 0). 
As the coupling is increased, the locations of the acoustic peaks are also shifted to smaller 
scales. However, the amplitudes of the odd-number peaks decrease due to the decrease of 
the baryon density at early time and the increasing ISW effect as the coupling scaled with 
the COBE normalization. The amplitudes of the even-number peaks increase due to the 
increase of the CDM energy density. Even though the ratios of the height of the first peak 
to those of higher peaks between different couplings are quite similar to one another, these 
are quite different from the ACDM model [13]. The location of the turnover in the matter 
power spectrum corresponds to the scale that entered the Hubble radius when the universe 
became matter-dominated (a eq ). Thus the coupling might shift the turnover scale. 

This paper is organized as follows. In the next section we show the basic equations of 
linear perturbations of the coupled quintessence model. In Sec. 3, we derive the formal 
entropy perturbation due to multiple fluids. We also consider the possibility of isocurvature 
perturbations of the quintessence. We check the coupling effects on CMB and matter power 
spectrum in Sec. 4. The effect of coupling on the metric perturbation is considered in Sec. 
5. Our conclusion is in the last section. 



2 Linear perturbations 

We will consider the perturbation effect of the quintessence model which is coupled to the 
CDM. First we start from the metric in the conformal Newtonian (longitudinal) gauge, 
which is restricted only for the scalar mode of the metric perturbations [14]. The line 
element is given by 



ds 2 = a 2 (rj) 



-(l + 2ty(ri,x))dr] 2 + (l - 2<$>{i 1 ,x))dx i dx l 



(2.1) 



where rj is the conformal time, \l/ is the amplitude of perturbation in the lapse function, 
and $ is the amplitude of perturbation of a unit spatial volume. We will consider only a 
flat universe case. If the coupling is derived by a Brans-Dicke Lagrangian, the radiation is 
decoupled from the dark energy [15]. Due to the strong constraint on the coupling to the 
baryons from the local gravity experiments such as radar time-delay measurements, we will 
assume that the baryons are decoupled from quintessence [16]. Also from the violation of 
the weak equivalence principle, even though it is in a still way that is locally unobservable, 
we can have the species-dependent couplings [17]. So we will make an assumption that 
the scalar field is coupled only to CDM by means of a general function exp[_B c (0)] and 
there is no coupling to the baryons or the radiation [18]. We can write the general equation 
including this interaction as 

S = -J d A x^{^-[d^d^ -R] + V{<t>) -C c -C r - C b \, (2.2) 



2 



where M = l/V8nG is the reduced Planck mass and £jS denote respectively CDM, radia- 
tion, and baryons. If we regard the matter as a gas of pointlike particles with masses m c 
and paths x"(t), then we can write C c as 

C c = -^5(x-x c (t))(-9^ c x^ 2 , (2.3) 

where m c = e Bc ^m* (m* being a bare mass of the CDM) [19, 20]. Each fluid element has 
an energy-momentum tensor T^ u where (3 includes all of the species. The total energy- 
momentum tensor is covariantly conserved, however the energy-momentum transfer between 
CDM and quintessence is written as 

= 0, (2.4) 

V M 7f 7)i , = 0, (2.5) 
= Q {d> , (2.6) 

where 7 denotes baryons or radiation, d denotes CDM or quintessence, and Q^)u is the 
energy-momentum transfer vector which shows the energy transfer to the (i-fluid [21]. This 
transfer vector is constrained as 

E^ = °- ( 2 - 7 ) 

d 

From the above action (2.2) we have the following equation which will give the constraint 
equation for the interaction between the quintessence and the CDM: 

From this we can find the unperturbed part of the field equation, 

+ m<p + w^f- = -w^r pc) { ] 

where the prime means d/drj, 7i = a' /a, and p c means the energy density of CDM. As we 
mentioned before the energy-momentum of each species may not be conserved due to the 
scalar field coupling even though the total energy momentum does conserve. If we use this 
fact, then we can rewrite Eqs. (2.4), (2.5), and (2.6) as 

Ptot = -3W(ptot+Ptot), (2-10) 

p 7 = -3H(p 7 +p 7 ), (2.11) 

p' c = -3n(p c + Pc ) + B c ^'p c = -3n(p c + p c )(l-B c ), (2.12) 

p's = -3H(p<t, +P4,) - B^4>'p c = -3H(p<t> - ( 2 - 13 ) 
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Figure 1: (a) Cosmological evolution of the equation of state parameter, u>^, and the energy- 
density parameters, f^, of each component for A = 5 and n c = 10~ 2 . The evolution of the 
parameters is similar for other choices of A. (b) Deviations of the scaling of the CDM energy 
density (£) from a -3 for couplings with n c = 1CT 3 (dashed line) and n c = 1CT 2 (dash-dotted 
line) . 



where B c ^ = dB c /d(p. Henceafter we will adopt the potential and the coupling as given in 
Ref. [20]: 



exp[B c (0)] 



Voexp 



A0 5 



V{<t>)/Vo X 



with b c + 1 > 0, 



V 1 + b c J 

where Vq, A, b c , and n c are constant parameters. Here we will consider only b c = case 



(2.14) 
(2.15) 



From these we have the evolution of the background quantities as shown in the first panel 
of Fig. 1. The parameters we use in this figure for the present energy density contrasts of the 
quintessence and the CDM are = 0.76 and = 0.20, respectively. To be compatible 
with observational data, the energy density of the quintessence must be subdominant during 
the big-bang nucleosynthesis (BBN) fi^ BN (a ~ lO" 10 ) < 0.2 at T ~ 1 MeV [10], and 
the energy density at present must be comparable to the preferred range for dark energy, 
0.60 < n^\a = 0) < 0.85 [2]. The stronger bound on the energy density of the quintessence 
during BBN, ^ BN (a ~ 10~ 10 ) < 0.045 [22], is also well satisfied in this model. We note 
that the equation of state (EOS) of the quintessence = p^/ ps u r = 1/3 during the 
radiation dominated era as in other tracker solutions. Instead of falling from the tracking 
value of 1/3 towards —1, as in the uncoupled case, ui^ increases towards higher values 
(< +1), before ultimately dropping to —1 at later times. This is due to the fact that when 
the matter density increases, the effect on the effective potential V e //(0) = V(<fi) + p c given 
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in Eq. (2.9) drives the faster evolution of 0. From Eq. (2.12), we find that the scaling of 
the CDM energy density differs from the usual a~ 3 due to the coupling : 

p c (a) = p^a- 3 (eMB c (m)-B c (m)])^P { c ) ^, (2.16) 
Z-x = B c (<l>(x))-B c (<i>(0)), (2-17) 

where we set the present scale factor as one (a^ = 1), x — In a, and £ is the deviation of the 
CDM redshift as a result of the coupling to a scalar field. This is shown in the second panel 
of Fig. 1. As the coupling is increased, the magnitude of £ is increased. The £ depends 
also on the scale factor. Consequently, the location of the turnover in the matter power 
spectrum will be shifted due to the coupling as well (see below). 

We can find the time-component of the source term in the unperturbed background 
(Q(d)o) from Eqs. (2.6), (2.7), (2.12), and (2.13) which is given by 

Q(c)o = -Q(ct>)o = -B c ^(p'p c . (2.18) 

Due to the conservation of the total energy- momentum (2.4), we can find the constraint 
equation of B c and B^, 

(p c + p c )B c +(p +p )i^ = O. (2.19) 
To consider the perturbation we can decompose the scalar field as 

<i>(ri,g) = $(ri) + 6<f>(ri,g), (2.20) 

where is the unperturbed part and 5<p is the perturbed part of the scalar field. We will 
express the perturbed parts of each quantities by means of Fourier expansions. So the above 
perturbed scalar field will be expressed as 

5<j ) ( V ,x) = J2SMv)e rk - £ , (2.21) 
k 

where 

SMV) =yj tyfa, x)e' ils d 3 x. (2.22) 

Therefore, the energy-momentum tensor of the scalar field can be decomposed into the 
unperturbed part and the perturbed one. The background energy-momentum tensors are 

M 2 - 

(0) T ( V = -(—fi 2 + V{4,)) = -p+, (2.23) 



^2a 2 

= (^> 2 -V(t))6i=p,6i = u;,p,6i, (2.24) 
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and the first-order perturbed parts are 

*T&o = ^(Mff-MW-« 2 ^), (2.25) 

^ = i(-M 2 ^* + M^V-a 2 ^^)^ (2.26) 

M 2 - 

= --^W). (2.27) 

We can repeat the similar consideration for the other components. We will regard the 
CDM as a perfect fluid of energy density p c and pressure p c . To the linear order in the 
perturbations the energy momentum tensors are given by 

T ( ° c)0 = -{pc + ^jp-8<l>Pc + Sp^j = - [pc + 5p)j , (2.28) 
T? c)l = (p c + p c )v (c)t = -T{ c)0 , (2.29) 

TUi = [pc + + *Pc)*j + % )3 = {pc + Sp^Si + Ej cy , (2.30) 

where v l = dx l /di] and E* is an anisotropic shear perturbation. Here we define the per- 
turbed part of the CDM as 5p c = B ^5<f>p c + 5p c , where the first term of the right hand side 
is due to the coupling of the scalar field to t the CDM 1 . For the baryons and the radiation 
we have the following equations, 

T(° 7 )o = -{Py + Spy), (2.31) 
= (p 7 +p 7 ) %) , = -T^ )0 , (2.32) 

T(y)j = (p T + 5p 7 )<5j + £j 7)j . (2.33) 

For a flat Friedmann-Robertson- Walker universe we can collect the unperturbed equations 
for each species: 

^ = ^{Pr + Pb + Pc + P<l>) = ^J2PP= J^Pcr, (2.34) 



M 2 V^' " " / t, r VJ - M2 z^rp- M2 , 



(2.35) 



p' 7 = -3H(p 7 +p 7 ), (2.36) 
p' c = -3H(p c +p c ) + B c ^'p c = -3H(pc + Pc)(l-B c ), (2.37) 
pi = -ZH{pt+Pt)-B^'p c =-m(p4+p4)(l-B4). (2.38) 



'We can show this as follows: p c = e B ^p* c = (e B ^ +B c {<j>) ^6<pe B ^)(p* c +5p* c ) ~ p c + B c (<P) ^6<pp c + 
8p Cl where we means B c (4>)^ = dB c {4>)/d<j> and p* is the bare energy density of the CDM. 



6 



If we include terms up to the first order of the energy transfer vector (2.6), then we can 
write Q( d ) v as 

Q(d) = -0(d)(l + *)-*Q(d), (2.39) 
Q(d)i = (f(d) + Q(d)V(d)i) .• (2.40) 

where 5Q(d) and f^) is the energy and momentum transfer of the CDM or quintessence, 
respectively. These terms should be included in the coupled quintessence models when we 
use the conformal Newtonian gauge [21]. If we missed this term, then we would have 2^B C ^ 
instead of 3^/B c ^ in the last term of Eq. (2.41). With using Eqs. (2.39) and (2.40), we can 
find the perturbed part of the energy momentum conservation equations in A;-space. The 
perturbed equations for the scalar field and the fluids which can be obtained from Eqs. (2.5) 
and (2.6) are given by 

Sf + ^ + 2^0' + 2*^Vj + ^V M 5cf>- (V + 3&)fi 

= -|U [ B cH^ + + + MB Ct ^ 

= ~^Pc (B cM 5cf> + B C /5 C + 3*B C ^ , (2.41) 

M = _(i + ^ c){ e c -W)-m( 5 -^ + l)5 c + B c -^{5 c + ^), (2.42) 

Pc OPc 

M = _ {l + u c ){d c -W)-m(^ + i)5 c + B cM $5<p 

Pc OPc 

+ B^05 c + 5<l>' + 4t*), (2.43) 
5' c = -(1+u c )(9 c -3&)-3H(^-u; c )5 c + B c ^, (2.44) 

K = ~(l+u; c )(9 c -3^)-3n(^-u; c )l + B cM m + B c ^ ( f ) ' 

+ B C> $'V, (2.45) 

9> c = -n(l-3u c )9 c --^9 c + S ^^k 2 5 c -k 2 a c + k 2 ^ 

1 + uj c 1 + uj c 

= -n(l-3cu c )e c -^^9 c + 6 -^l^k 2 5 c -k 2 a c + k 2 ^ 

1 + UJ C 1 + tu c 

' B.$'8 C , (2.47) 



where we define S c = 5p c / p c , 5 C = 5p c / p c , 9 C = ik ■ v, cu c denotes the equation of state 
parameter (EOS) of the CDM, and a c is related to the CDM anisotropic stress perturbation 
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II C , by <7 C = 2n c p c /3(p c +p c )- We express both ^ and S' c explicitly in order to show that 
the equations of the energy density perturbation can be different with different definitions. 
However, it is the coupled energy density which is measured in observations. As such, we 
will use 5p as the energy density contrast of each species. Every term containing B c (<f>) 
in the above equations comes from the coupling. If we drop all these terms, then they 
are obviously identical to the expressions Ref. [23]. If we have more than one ideal gas 
components, then we have entropy perturbation (SS) which can be expressed as 

tiptat = ^2[(dpp/dpp)\sSfhat + (dpp/dSf3)\ p 5S tot = c 2 tot 5p tot +Ptotr to t, (2.48) 

P 

where c 2 ot is the overall adiabatic sound speed squared and r to t is the total entropy per- 
turbation. We will consider this more carefully in the following section. We can write the 
perturbed equations for the CDM from the above generic perturbation equations (2.45) and 
(2.47). However due to the coupling between the baryons and the photons we also need to 
consider the Thomson scattering term in the baryon-photon fluid. After including this we 
have the following equations. 

5' b = -6 b + 3&, (2.49) 

9' b = -H9 b + c 2 k 2 5 b + k 2 ^ + ^an e a T (9 1 -9 b ), (2.50) 

S' c = -6 c + 3^ + B^'% (2.51) 
S' c = -9 C + W + B cM 4>'5<t> + B c j5<t>' + B C ^, (2.52) 
0' c = -<H9 c + k 2 V, (2.53) 

K- = ~9 r + W, (2.54) 

9' r = k 2 (^5 r - a r ) + k 2 ^ + an e a T (9 b -9 1 ), (2.55) 
where n e is the electron number density, ot is the cross section for the Thomson scattering. 



3 Entropy perturbation 

Let us start from the definition of the total energy density perturbation and that of the 
total pressure perturbation, 

tiptot = ^Spp. (3.2) 

P 

For a given ptot(p, S), the pressure fluctuation can be expressed as 

Mot = Ctot^ptot + ptotr in t + ptotr re i, (3.3) 
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where S is an entropy, c 2 ot is the overall adiabatic sound speed squared, r int and r re i are 
the intrinsic and the relative entropy perturbations respectively. We have 

Ptotr int = ^ppTp, (3.4) 



Pto t r rc i = ^2(c} - c 2 tot )5p/3, (3.5) 

P 

4, = (3-6) 

Ptot 

where T^t is the sum of the intrinsic entropy perturbation of each fluid and r re i arises from 
the relative evolution between fluids with different sound speeds. As we mentioned before 
the energy momentum of each species may not be conserved due to the scalar field coupling 
even though the total energy momentum does conserve. We can rewrite the equation for 
the total adiabatic sound speed (3.6) as 

c t ot = 2^ W - B P) ~ —Z • ( 3 - 7 ) 

Ptot + Ptot 

Now we can rewrite the relative entropy perturbation as 

Pt ot r rcl = \ _ c l )Spa + £ Bpc% P p + P p)A, (3.8) 

1 p, a Aot + Ptot p 

where 

Sp a = Ap - A a , (3.9) 
A, = (3.10) 

P/3+P/3 

A = 5ptot , (3.11) 

Ptot + Ptot 

where Sp a is the entropy perturbation [21]. Due to the A-term the relative entropy pertur- 
bation is non-vanishing even without the non-adiabatic perturbation. This is improper, so 
we redefine the new quantities as 



Ae = (i-JWto+w)' (3 ' 12) 

V = A„-A,,. (3.13) 
With these quantities we can rewrite Eq. (3.8) as 

„ r _ 1 \- (I - Bp)(l - B a )( P p + P p)(p a + Pa ) 2 2 ~ 

PtotJ- rel - - 2^ , l C /3 ~ C a) b Pa- (A i4 J 

z p jCt Ptot T Ptot 
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Figure 2: (a) Cosmological evolution of the adiabatic sound speed c| of the scalar field for 
different values of n c = (dotted line), n c = 1CT 3 (dashed line), and n c = 1CT 2 (dash-dotted 
line) when A = 5. (b) Behavior of the coupling modification, B^ as a function of the scale 



factor a for different values of n c . 



3.1 Isocurvature condition 



Due to the out of thermal equilibrium nature of the quintessence, we need to check the 
isocurvature evolution of the scalar field perturbation [24]. We can analytically show this 
in the tracking region. First of all the adiabatic sound speed squared c| of the quintessence 
can be represented as 



1 + 



3W(1 + ^)(1-B^) 



(3.15) 



where we have used Eq. (2.13). From this equation we can find the second derivative of the 
potential, 



M 2 



L^-i-|(4+ 1 )( 1 -^)]-f w (4- 1 ) fl i- 

(3.16) 

The relation between Ti, 2 and Ji! can be found from Eq. (2.35) : 



T~C' 1 

y2 = -^(^^Ot), 



^2 U 13^/3^ 



(3.17) 
(3.18) 



where u tot is the weighted EOS. As long as we use the background as shown in the pervious 
model, the tracking region is well established during the radiation dominated era [20]. 
During this era we can find the following relation, 



iO^ — u r . 



(3.19) 
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This is shown in the first panel of Fig. 2. As the coupling is increased we can have the 
non-monotonic behavior of as shown in the n c = 1CT 2 case. In addition, we can rewrite 
the coupling term as 

l>,,,o'l>, H 



Mto+p*) 4>~ c '«"" c - (3,20) 

As shown in the second panel of Fig. 2, this term is negligible during the radiation dominated 
epoch. The coupling drives a faster evolution of when matter energy dominates the 
Universe and the magnitude of 3$ depends on the energy density of the CDM. With these 
facts we can have the approximate expression of Eq. (3.16), 



M 2 



V: 



H 2 (d 



n 2 



3cu tot + 3cl + 4 



(3.21) 



where we have used Eq. (3.17) in the second equality. The isocurvature mode in the 
radiation dominated epoch can be obtained from Eq. (2.41) and using the fact that a oc rj 
during the radiation dominated era. To check this we can put $ = and then we have 



+ 



(3ccv + l)r] 



+ 



k 



1) 



(3oo r + l) 2 r,' 



-(3u; r + 3cJ + 4) 



0. 



If we use the tracking solution (3.19), then we can rewrite the above equation as 

= 0. 



V 



+ 



,2 3 
k 2 + — 



(3.22) 



(3.23) 



The solutions of this equation are Bessel functions, 

50(77) = const ?T 1/2 J±\iVTi/2\(kri)- (3-24) 

Thus both solutions decay in time. At the superhorizon scale the /c-dependent term can 
be neglected and we obtain the power-law solutions, 5 (ft oc i] u , where the power index v = 
(— 1 ± i^/TT)/2. Therefore, any initial nonzero isocurvature fluctuation of the quintessence 
is damped to zero with time. We will use only the adiabatic perturbations in the following 
section. 



4 Effects of coupling 

The non-minimally coupled quintessence models have been investigated as a possible solu- 
tion for the late time coincidence problem [11]. The coupling gives rise to the additional 
mass and source terms of the evolution equations for CDM and scalar field perturbations. 
This also affects the perturbation of radiation indirectly through the background bulk H 
and the metric perturbations [25]. The value of the energy density contrast of the CDM 
(fl c ) is increased in the past when the coupling is increased. 
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Figure 3: (a) CMB large-scale anisotropy power spectra of ACDM (solid line), minimally 
coupled n c = (dotted line), and non-minimally coupled n c = 10~ 3 , 1CT 2 (dashed, dash- 
dotted line respectively) quintessence models, (b) Same spectra for the entire scales. 



4.1 CMB 



The temperature anisotropy measured in a given direction of the sky can be expanded in 
spherical harmonics as 

AT 

© = -=-(n) = 22a em Ye m (n), (4.1) 
1 e m 

where is the temperature brightness function that is the fractional perturbation of the 
temperature of the photons T = T (l + 0), n is the direction of the photon momentum, 
and a£ m are the multipoles. We can also expand the brightness function as a Legendre 
polynomial, Pi : 

e = X)H) / (2^+1)6^/. (4.2) 

i 

From the inflationary scenario, the multipoles are Gaussian random variables which satisfy 

{a* im ai> m >) = Ci 5 U i 8 mm >. (4.3) 

The angular power spectrum (CV) contains all the information about the statistical prop- 
erties of CMB and is defined as [26] 



4tt 2vr 2 J 1 k 2£ + 1 

In the standard recombination model the acoustic oscillations will be frozen into the 
CMB. A generalization of the free-streaming equation in a flat universe gives the resulting 
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anisotropies 



[0 O + (msMHv - vis)) + QiM^MHv - m*)) + J" (*' - $').MM?7 - *jQ)*7 

(4.5) 



2£+ 1 



where j£ is the spherical Bessel function and r]i s is the conformal time at last scattering. 
Photon density perturbation is related to the temperature perturbation in the matter rest 
frame : 



and the gravitational potentials will be given in Eq. (5.1) - (5.4) in the following section. 
We use the fact that in the absence of anisotropic stress, the two scalar potentials \& and $ 
defined in the conformal Newtonian gauge (2.1) are equal and they coincide with the usual 
gravitational potential in the Newtonian limit. 

Now, we investigate the effects of non-minimal coupling of a scalar field to the CDM 
on the CMB power spectrum. Firstly, the Newtonian potential at late times changes more 
rapidly as the coupling increases, as shown in Eq. (5.1). This leads to an enhanced ISW 
effect as indicated in the last term of Eq. (4.5). Thus we have a relatively larger Ce at 
large scales (i.e. small £). Thus, if the CMB power spectrum normalized by COBE, then 
we will have smaller quadrupole [27]. This is shown in the first panel of Fig. 3. One 
thing that should be emphasized is that we use different parameters for the ACDM and 
the coupled quintessence models to match the amplitude of the first CMB anisotropy peak. 
The parameter used for the quintessence model is indicated in Fig. 3 (i.e. = 0.76, 

Q$ = 0.191, = 0.049, and h = 0.7, where h is the present Hubble parameter in the 
unit of 100kms _1 Mpc _1 ). However, these parameters are well inside the 1 a region given by 
the WMAP data. We use the WMAP parameters for the ACDM model (i.e. = 0.73, 
= 0.23, fi[ 0) = 0.04, and h = 0.72). In both models we use the same spectral index 
n s = 1. The heights of the acoustic peaks at small scales (i.e large £) can be affected 
by the following two factors. One is the fact that the scaling of the CDM energy density 
deviates from that of the baryon energy density. Therefore for the given CDM and baryon 
energy densities today, the energy density contrast of baryons at decoupling (fi^) is getting 
lower as the coupling is being increased. This suppresses the amplitude of compressional 
(odd number) peaks while enhancing rarefaction (even number) peaks. The other is that 
for models normalized by COBE, which approximately fixes the spectrum at £ ~ 10, the 
angular amplitude at small scales is suppressed in the coupled quintessence. This is shown 
in the second panel of Fig. 3. The third peak in this model is smaller than that in the 
ACDM model. The WMAP data do not show the value of the third peak but quote a 
compilation of other experiments [28]. The ratio of the amplitude between the second and 
the third peaks is 1.03±0.02. In the ACDM model this value is 0.986 and in our model these 
values are 1.08 and 1.11 for without and with the coupling equal to n c = 0.01 respectively. 



5 r = 4© + m^, 



(4.6) 
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We also show the n c = 0.02 case in this figure. With the same parameters this case can be 
ruled out from the current data. In all cases, we have used the CMBFAST code [29] with 
the modified Boltzmann equations to compute the CMB power spectrum. 

Secondly, for £ > 200 we can see that the locations of the acoustic peaks are slightly 
shifted to smaller scales (i.e. larger £). This can be explained as follows. The locations of 
peaks and troughs can be parametrized as 

t m = £A(m-(p m )=£ A (m-(p-6(p m ), (4.7) 

where £a is the acoustic scale dependent on the geometry of the Universe, tp is the overall 
peak shift, and 5<p m is the relative shift of the m-th peak relative to the first one [6]. The 
overall peak shift, <p is given by 

^ .267(^) ai , (4.8) 

where r is = p^/pffl is the ratio of the energy densities of radiation to matter at last 
scattering. The shift is due to both driving effects from the decay of the gravitational 
potential and contributions from the Doppler shift of the oscillating fluid. The acoustic 
scale £a depends on both the sound horizon si s at decoupling and the angular diameter 
distance D to the last scattering surface: 

D rj - rj ls 

£a = 7T = 7T— , (4.9) 

sis c s r]i s 

where c s is the average sound speed before last scattering : 

c s = rjZ 1 r U c s dr] with c7 2 = 3 + -— . (4.10) 
Jo 4 p r 

Also from the Hubble parameter, 

=H 2 Jn^ + n^a + n^a^ + ^-aA, (4.11) 



drj 



joy 

fJcr 



where Hq is the present value of the Hubble parameter, we can find the angular diameter 
distance, 

r a<-°) A n 

D = Vo -Vis = H 1 . (4.12) 

Ja ls ^(0) + n (0) a + Q (0) a l +e + p / p (°) a 4 

The sound horizon is not affected by the coupling and the effect of coupling on the angular 
diameter distance is also quite small in our model. However, the overall shift and the relative 
shift are affected by the coupling. As the coupling is increased, the value of r is is decreased. 
Hence, the locations of the peaks are shifted to the right. However, this shift is quite small. 
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Figure 4: Matter power spectra for the models using the same parameters in Fig. 3. 

We show this in the second panel of Fig. 3. The heights and the locations of the acoustic 
peaks in various models are summarized in Table 1. There is a significant difference for the 
heights of the peaks of the second and the third peaks between the models. Thus upcoming 
observations continuing to focus on resolving the higher peaks may constrain the strength 
of the coupling. 

Table 1: Summary of the heights and the locations of the considered models. The heights 
of peaks, Ai have been scaled by a factor of 10 10 . A i: j is the ratio of the i-th peak height 
to the j-th. The locations of peaks, £{. 



Model 


A, 


A 2 


^3 


ii 




h 


A V .2 


^2:3 


ACDM 


7.26 


3.25 


3.30 


220 


536 


812 


2.23 


0.986 


n c = 


7.04 


3.23 


3.00 


217 


532 


811 


2.18 


1.08 


n c = 1Q- 2 


7.02 


3.36 


3.03 


218 


535 


818 


2.09 


1.11 



4.2 Matter power spectrum 



The effects of quintessence on structure formation in several other models have been inves- 
tigated [30, 31]. Structure as a function of physical scale size is usually described in terms 
of a power spectrum : 

P(k) = {\5 k \ 2 )=Ak n *T 2 (k), (4.13) 
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where A is the COBE normalization, n s is a power index, and T(k) is the transfer function. 
The coupling of quintessence to the CDM can change the shape of matter power spectrum 
because the location of the turnover corresponds to the scale that entered the Hubble 
radius when the Universe became matter-dominated. This shift on the scale of matter and 
radiation equality is indicated in the second panel of Fig. 1: 

a eq ~ ^jexpb3 c (0 o ) - B c (<f> eq )}, (4.14) 

Pc 

where pf> and p^ are the present values of the energy densities of radiation and CDM 
respectively, and the approximation comes from the fact that the present energy density of 
CDM is bigger than that of baryons (p^ > pf^). Increasing the coupling shifts the epoch 
of matter-radiation equality further from the present, thereby moving the turnover in the 
power spectrum to smaller scale. If we define k eq as the wavenumber of the mode which 
enters the horizon at radiation-matter equality, then we will obtain 

k eq = — . (4.15) 

Veq 

However, from the previous subsection, we notice that the value of i] eq remains unchanged 
for different couplings and this degeneracy is indicated in Fig. 4. We have used different 
parameters for the ACDM and quintessence models; the matter power spectra look different 
between the models. There is a slight suppression in the quintessence models. Note that 
a bias factor could resolve the discrepancy and perhaps a parameter fitting may also help 
this. However, the detailed parameter fitting is beyond the scope of this paper. 

We can write the equation of the matter fluctuation in the synchronous gauge during 
the matter dominated epoch : 

K + m> - fo Ofc + fo) _ (^W = o, (4.i6) 

2 p cr a 



where the coupling term F((j>) is given by 

F(<f>) = B cM 4>'5<t> + B^54>' = n c A(0<50)'. (4.17) 
This equation can be rewritten by the structure growth exponent /, which is defined as 

Hence, we have the following equation for f(a) which is identical to Eq. (4.16): 
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If we use Eq. (2.16), then we can find that a oc rj 2 ^ 1 '^ and H — 2/(1 — £)?? 1 f° r a matter 
dominated universe. As we can see from the above Eq. (4.19) that the coupling term is 
negligible because varies much slower than Ti. Thus we can ignore the last term in 
Eq. (4.16). From this we can rewrite the above equation in a matter dominated era : 

5 " + (i4)¥ 4 = (4 ' 20) 

This equation has two solutions, 

5f = c ±V ^, (4.21) 

where c± are arbitrary constants and 



- i + 0± V 24 + 1 + 2 , x 

v± = — — — + r— : — . 4.22 

2(1-0 1 ' 

5+ indicates a growing mode, which is only relevant today because |<5+| is small at early 
time and we can ignore the decaying mode 8~. If we remove the coupling effect in this 
solution, then we can recover the well known solution 5+ = rf^. This effect is shown in 
Fig. 4. 

As the coupling is increased, we have little more matter at early times and this increases 
the height of the matter power spectrum. Again this effect is tiny and we can hardly see 
the difference between various couplings. 



5 Metric perturbation 

In addition to these, the perturbed equations of the metric can be obtained from the Einstein 
equations: 

k 2 <$> + m(& + mi) = %— 5T ( L)o, (5.1) 

^ Ptot 

k 2 1 Ti 2 

<f>" + W + 2$') + (27i! + H 2 )^ + —(<$>-■$>) = -—5T{ tot)i , (5.2) 

k 2 U' + ny) = i^W 2 £(1 + ^)^3, (5-3) 

k 2 ($-y) = ^£(1+^)^, (5.4) 

where we have used the unperturbed equations (2.34) and (2.35) . From Eqs. (5.1) and 
(5.2), we can find the metric perturbation equation, 

k 2 

$" + H(^' + 5$') + 2(H' + 2H 2 )^ + —(4$ - 

3 

a 2 

[(1 - C tot)$Ptot ~ PtotTint - PtotTrel]- (5.5) 



2M 2 
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If we use Eq. (5.1), this equation can be rewritten as 



v' + (2 + 3cL)<f'J + [m 1 + n\\ + 3cL) 

2 

2M 2 



3 



l+34j$-tt 

2 



— 7}~j\^2 (PtotF int +Ptot^rel S ) = 2 Ptot^tot, (5-6) 



where we define r tot = Y int + T re ; in the last equality. The last term in this equation comes 
from the coupling of the scalar field. So even if we start from the adiabatic condition ( 
Ptot^tot = 0), we can analytically solve the above equation for the specific case. Let us first 
put ^ = $ (no anisotropic stress) and k 2 c 2 ot § = (consider the superhorizon scale), then 
Eq. (5.6) is simplified as 



.2 N 



2M 2 



Ptot^tot- 



$" + m{\ + c 2 j$' + [2H' + n 2 (i + 3cf 

If you use the background equations, this equation can be rewritten as 

m 2 Ptot ] 

2 ptot 



$" - (\n[p tot +ptot])'^ + (HPtot/(ptot + Ptot)])' HQ 



tot- 



(5.7) 



(5.8) 



If we consider the adiabatic condition (ptot^tot — 0), then the above equation is identical to 
the equation in Bardeen's article [32]. As such, the metric perturbation can be rewritten as 
the curvature perturbation [33], 



(H-^' + Q) 3^ 
(l + wtot) 2 

With this we can express Eq. (5.7) in the adiabatic case as 



(5.9) 



C = o. 



(5.10) 



This result looks the same to the minimally coupled case and there seems to have no 
difference from the non-minimally coupled case. Nevertheless, the coupling information is 
absorbed in both c 2 ot and the perturbation equation of each species. 



6 Conclusions 

We have analyzed the linear perturbations of the cosmological models for the scalar field 
with its self-interaction potential, V((f>) = Vo exp(A0 2 /2), and its coupling to the CDM, 
exp[_B c (0)] = exp(n c A0 2 /2). The evolution of the non-perturbed background scalar field 
4> occurred in the tracking regime throughout the radiation dominated epoch. The full 
analysis of the non-perturbed and perturbed equation of each species in the conformal 
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Newtonian gauge has been done including the proper energy-momentum transfer vector 
due to the coupling, Q(d)w The Boltzmann equations have been modified as to account for 
the coupling between a scalar field and the CDM. 

We have seen that the energy-momentum of each species may not be conserved as a 
result of the coupling of the scalar field to the CDM. Thus we have redefined the con- 
cepts of entropy perturbations. We have shown that the isocurvature perturbation of the 
quintessence has been damped to zero with time in the tracking regime during the radiation 
epoch. Thus we have constrained our considerations on the adiabatic perturbation. 

We have considered the CMB anisotropy spectrum and the matter power spectrum for 
the non-minimally coupled models. Additional mass and source terms in the Boltzmann 
equations induced by the coupling give the rapid changes of the Newtonian potential $ and 
enhance the ISW effect in the CMB power spectrum. The modification of the evolution of 
the CDM, p c = p^a~ 3+ ^, changes the energy density contrast of the CDM at early epoch. 
We have adopted the current cosmological parameters measured by WMAP within la level. 
With the COBE normalization and the WMAP data we have found the constraint of the 
coupling n c < 0.01. The locations and the heights of the CMB anisotropy peaks have been 
changed due to the coupling. Especially, there is a significant difference for the heights of 
the second and the third peaks among the models. Thus upcoming observations continuing 
to focus on resolving the higher peaks may constrain the strength of the coupling. The 
suppression of the amplitudes of the matter power spectra could be lifted by a bias factor. 
However, a detailed fitting is beyond the scope of this paper. The turnover scale of the 
matter power spectrum may be used to constrain the strength of the coupling n c . 

Finally, we have investigated the metric perturbations including the coupling between 
the scalar field and the CDM. There is no difference to the curvature perturbation ( for 
the different couplings in the adiabatic case. However, the effects of the coupling have been 
absorbed in the Boltzmann equations already. 
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